Small Shear Viscosity of a Quark-Gluon Plasma Implies Strong Jet Quenching 
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We derive an expression relating the transport parameter q and the shear viscosity 77 of a weakly 
coupled quark-gluon plasma. A deviation from this relation can be regarded as a quantitative 
measure of "strong coupling" of the medium. The ratio /q, where T is the temperature, is a 
more broadly valid measure of the coupling strength of the medium than 77/s, where s denotes 
the entropy density. Different estimates of q derived from existing RHIC data are shown to imply 
radically different structures of the produced matter. 
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The highly excited, strongly interacting matter formed 
in collisions of large nuclei at the Relativistic Heavy Ion 
Collider (RHIC) exhibits two miusual properties: The 
emission of hadrons with large transverse momentum is 
strongly suppressed in central collisions [H and the col- 
lective flow of the matter is well described by relativistic 
hydrodynamics with a negligible shear viscosity [2| . The 
suppression of hadrons at large transverse momentum, 
generally referred to as jet quenching, is understood to 
be caused by gluon radiation induced by multiple col- 
lisions of the leading parton with color charges in the 
near-thermal medium 0, 0, @] ■ This process is governed 
by a transport parameter q, which is the squared aver- 
age transverse momentum exchange between the medium 
and the fast parton per unit path length [3| ■ 

The RHIC data have been interpreted to imply that 
the quark-gluon plasma produced in the nuclear collision 
is a strongly coupled medium, which may not even con- 
tain quasi-particles whose interactions can be treated in 
some effective perturbation theory Q. While the abso- 
lute values of the shear viscosity-to-entropy ratio 77/5 de- 
termined from model analysis of RHIC data tend to rule 
out a class of weak coupling approaches based on the 
Hard-Thermal Loop (HTL) approximation, it may still 
be consistent with the general picture based on weakly 
coupled quasi-particles with partonic (quark and gluon) 
quantum numbers. The present letter outlines a method 
for quantitatively testing this hypothesis. 

The interaction of hard jets with strongly interacting 
matter may always be treated in a perturbative expan- 
sion involving short lived partonic excitations as the QCD 
coupling is weak at short distances. Thermal excitations 
in a QCD medium, on the other hand, may or may not be 
partonic quasi-particles. Here, we take the term "weakly 
coupled" to mean that the properties of the medium 
can be described pcrturbativcly on the basis of an ap- 
propriate, partonic quasi-particle picture. This notion 
does not preclude the possibility that the quasi-particles 
themselves may emerge nonperturbatively from the fun- 
damental quanta of QCD. "Strong coupling" indicates 
the absence of such a partonic quasi-particle description. 
The central assertion advanced here is that there exists 



a general relation between the jet parameter q and the 
coefficient of shear viscosity 77 which holds for any weakly 
coupled partonic plasma where the interaction between 
the quasiparticles, which is responsible for the genera- 
tion of a viscosity, has the same structure and strength 
as the interaction between the leading jet parton and the 
medium. 

In any partonic quasi-particle framework, the trans- 
port parameter governing the radiative energy loss of a 
propagating parton in SU(3)-color representation i? is Q : 
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The shear viscosity of a fluid is defined as the coeffi- 
cient of the contribution to its stress tensor, which is 
proportional to the divergence-free part of the velocity 
gradient. In the framework of kinetic theory based on a 
quasi-particle picture, the shear viscosity 77 is determined 
by the mean-free path At(p) of a constituent particle of 
momentum p in the medium: 



77 ~ Cp{p)Xf, 



(2) 



with C « 1/3 0, i]. A heuristic connection between 
77 and q can be established by the observation that the 
mean-free path is related to the average transport cross 
section of a quasi-particle in the medium: Af = (pcrti )^^. 
When soft scattering dominates, as in the case of pertur- 
bative QCD, the transport cross section is related to the 
differential cross section by the relation: 
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where is the center-of-mass energy. For a thermal 
ensemble of massless particles, (p) « 3T and (s) « 18T^, 
and thus: 



77 ; 
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Using the relation s ~ 3.6 p for the entropy density of a 
gas of free, massless bosons, the following heuristic equa- 
tion relating the ratio of the shear viscosity to the entropy 
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density with the transport coefficient may be derived: 
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(5) 



This relation shows that a large value of q implies a small 
value for the ratio 77/s, which is thought to be bounded 
by the quantum limit 77/s > (47r)~^ [9|. 

When perturbation theory can be applied to the 
medium, the jet quenching parameter can be expressed 
in terms of the light-cone gluon correlation function 



^ dy- {Q)F^+{y'))e^P y\ (6) 



In Eq. dl]), (O) = (27r)-3 J d^p/2p+ f{p){p\0\p) denotes 
the ensemble average of an operator O in the medium 
composed of states \p) with occupation probability f{p), 
^ = (fc|n) /2£'(p+), (fc|n) is the average transverse momen- 
tum carried by the gluons in \p), and the index i sums 
over the two transverse directions, p = J d^ p f {p) / {2t:Y 
denotes the density of scattering centers, mainly gluons, 
in the matter. dafi/dq\ is the differential cross section 
for a parton on a scattering center. A nonperturbativc 
definition of qn in terms of a Wilson loop along the light- 
cone has been proposed in ref. • 

In principle, qR{E) depends weakly on the energy E 
and virtuality of the jet-inducing parton. In the limit of 
a thick (large opacity) medium of thickness L, the virtu- 
ality of the parton is determined by the saturation scale 
Qs of the medium , and the transport parameter for 
parton energy loss takes on the universal form q = Q^/ L. 
We will neglect a possible logarithmic energy dependence 
of q here. For a massless parton in a weakly coupled gluon 
plasma: 



qR 
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(7) 



in the leading log approximation ll|. In the remainder 
of this paper, we will omit the subscript in the notation 
for the transport parameter q = qA for a gluon jet. 

The expression for the shear viscosity depends on 
the mechanism that limits momentum transport in the 
medium. We start with the shear viscosity due to parton 
collisions in the quark-gluon plasma, following Arnold et 



al. |12l |. In the Chapman-Enskog approach to linearized 



transport theory, one parametrizes the deviation /i (p) of 
the parton distribution from equilibrium due by a func- 
tion A{p) in the form (using the notation of [13[) 



flip) = --^^P^P]{^u),j, 



(8) 



where {yu)ij denotes the traceless velocity gradient. One 
can then derive an integral equation for A(p) from the 
linearized Boltzmann equation. An analytic estimate 
can be obtained by restricting A to the functional form 



A(p) = Ap/T. The shear viscosity is expressed in terms 
of A(p) as 



15 T 



(9) 



where fo{p) denotes the equilibrium distribution of par- 
tons. Assuming that collisions in the medium arc domi- 
nated by soft scattering, the integral over the differential 
scattering cross section can be extracted from the colli- 
sion integral, yielding a factor q/ p, where q is defined by 
Eq. Il]) for a gluon jet. After a lengthy calculation one 
finds that in this limit the shear viscosity for a thermal 
gluon plasma can be expressed in the form: 



,T^p 



(10) 



In order to extract the constant, we can make use of the 
known result for the shear viscosity of a pure gluon gas 
in the leading logarithmic approximation [l^ : 



'7LL 



3.81 - 1 
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Inserting the jet quenching parameter in the weak cou- 
pling limit ([7]) and using the expression for the density 
of a thermal gas of free gluons, we obtain C" « 4.85. 
Comparing with Eq. (jj]), we find that this corresponds 
to C « 0.36. 

We next consider the case of anomalous shear viscos- 
ity, which is generated by dynamically created turbulent 
color fields in a rapidly expanding quark-gluon plasma 
[ist . Unstable collective plasma modes 14 1 result in the 
growth of long wave length modes of the glue fields, which 
ultimately saturate due to their nonlinear self-interaction 
[lil . [l6| . The turbulent quark-gluon plasma is character- 
ized by a non-vanishing expectation value of the gluon 
correlator {F'^+ {0)F+{y-)) of the soft color fields. 

The presence of localized domains of color fields in- 
duces an anomalous contribution to the shear viscosity 
of the matter 13[ . This contribution is derived by consid- 



ering the effect of random color fields on the propagation 
of quasi-thermal partons, which is described by a Fokker- 
Planck equation of the form: 



d_ 

dt 



V • - VpD^(p,i)Vf 



/(p,r,i) = (12) 



with the average parton phase space distribution / and 
the diffusion tensor 

Df,^ f dt' {FMt'),mim.t)) . (13) 

Here F = t/Q^ (E° + v x B°) is the color Lorentz force 
generated by the turbulent color fields and is color 
charge of the parton in the given representation. The 
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anomalous shear viscosity implied by Eq. (|12p can be 
evaluated, e. g., for randomly distributed color fields with 
a correlation time along the light-cone (isj : 



16C(6)(iV,^-l)^T6 



(14) 



The diffusion coefficient Dij in Eq. p^ . is closely re- 
lated to the transport parameter for radiative energy loss, 
Eq. Applying the definition of qa as rate of growth 
of the transverse momentum fluctuations of a fast gluon 
to an ensemble of turbulent color fields, one finds (for 
partons in the adjoint color representation): 



3(iV2 



1 



(15) 



Combining this expression with Eq. (|14p we obtain a re- 
lation between the contributions of turbulent color fields 
to the transport parameter in gluon energy loss and the 
anomalous shear viscosity of gluons: 



32C(6)(iV2 - 1)T6 



37r2 



4.51 



(16) 



The numerical coefficient here is very close to the one ob- 
tained for the coUisional viscosity. Alternatively, we can 
calculate the dimensionless ratio of the shear viscosity 
and the entropy density. For a free thermal gluon gas, 
the entropy density is given by 



s = 2{Nl - 1; 



45 



rp6 



(17) 



s 



which implies the following result for the shear viscosity- 
to-entropy ratio: 

— — «1.25— , (18) 
bi q q 

again in agreement with ^ for C ~ 1/3. This confirms 
our assertion that Eq. ([5]) holds if the transport prop- 
erties of a quark-gluon plasma are described by kinetic 
theory based on partonic quasi-particles. 

In order to explore the opposite limit, we first consider 
a solvable strongly coupled theory, = 4 supersymmet- 
ric Yang-Mills (SYM) theory. Analytical results for i]/s 
and q have been derived for this theory in the limit of 
large 't Hooft coupling A = g^Nc- The strong coupling 
expression for the ratio of shear viscosity to entropy den- 
sity is well established [isl . [lit : 



1 
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135 C(3) 
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(19) 



There is not yet universal agreement about the correct 
generalization of the jet quenching parameter. Here we 
follow Liu et al. [loj who defined q by means of an adjoint 
Wilson loop along the light-cone and found the result 



. 7r3/^r(3/4) /- 3 
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7.53 VAT^ 



(20) 



Inserting this expression into the right-hand side of ([5]) 
for C w 1/3 we obtain 



1.25- 



T3 



0.166 T] 

VA s 



(21) 



for A > 4.35. Obviously, the left-hand side of (|2T]) is 
a better measure of the strength of the coupling in the 
strong coupling limit than the right-hand side. 

The strong coupling limit of QCD exhibits quark con- 
finement. As is well known, this limit of QCD is also 
amenable to rigorous calculations in the chiral limit and 
at low temperature. In this limit, the shear viscosity- 
to-entropy density ratio of a pion gas is given by [22| : 



167rr4 



(22) 



The ratio grows rapidly with decreasing temperature, be- 
cause the pion-pion interaction weakens at low momen- 
tum. On the other hand, the jet quenching parameter 
for a gluon jet propagating through a pion gas has the 
form 



M _ 47r N,as 



(23) 



where G^(a;) is the gluon distribution function of the 
pion. Since the pion density in the chiral limit is pro- 
portional to T"^, this implies that the left-hand side of 
(|2ip is approximately independent of T: 



1.25- 



0.235 



9a 



x-tO 



(24) 



Again, /q is a much better measure of the underlying 
QCD coupling strength in this domain than t]/s. We note 
that T^/q <C rj/s in both cases of strong coupling. It is 
not surprising that the heuristic relation ([4]) does not hold 
in the strong coupling limit: In a strongly coupled sys- 
tem, thermal modes are not approximately described as 
elementary quasi-particles [2^ probed by jet quenching. 
On the other hand, highly energetic excitations can re- 
tain a quasi-particle nature, and their interactions with 
the medium continue to track the strength of the cou- 
pling. The violation of the relation ^ may thus be con- 
sidered as a general criterion for "strongly coupling" in 
a gauge theory. 

We can now ask whether the relation (j4|) applies to 
the matter produced at RHIC, according to the phe- 
nomenological analysis of jet quenching measurements. 
At present, the values of q deduced from the data by 
means of different approaches to jet quenching differ from 
each other by substantial factors [30|. Recent studies 
within the framework of the twist expansion, fitting to 
the experimental data on hadron suppression in the most 
central Au-fAu collisions at RHIC jisl . [26j gave values 
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1 — 2 GeV^/fm for the gluon quenching parameter 
at the initial time tq ~ 1 fm/c. A study using a variant 
of the eikonal approach [1] resulted in the estimate [l^l 
^o"^"^ = 10 — 30 GeV^/fm for the same quantity (scaled 
from the value given in p?} for a propagating quark). 
As we now show, the two different values imply radically 
different properties of the medium. 

One can estimate 28| the average initial entropy den- 
sity from the measured charge hadron multiplicity of 
the most central Au+Au collisions as sq = (33 ± 3) 
fm""^ at To = 1 fm/c corresponding to an initial tem- 
perature To = (337 ± 10) MeV. Using Eq. © with 
C = 1/3 one then finds that the lower value Qq^^^ im- 
plies 1.25 Tg /go = 0.12 — 0.24. whereas the higher value 
g^"*"^ implies 1.25roVgo = 0.008-0.024. This first result 
is close to the conjectured lower bound for rj/s and lies 
well within the range of values for the shear viscosity-to- 
entropy density ratio {rj/s < 0.3) which are compatible 
with the measured hadron spectra from Au-|-Au collisions 
at RHIC H^. If ultimately confirmed, it would indicate 
that the quark-gluon plasma produced in Au-|-Au colli- 
sions is (marginally) weakly coupled. On the other hand, 
the second result lies far below the lower bound for rj/s 
and thus implies that the quark-gluon plasma formed at 
RHIC is deep in the strong coupling regime and cannot 
be described as a quasi-particlc plasma. 

In summary, we have derived a general relation be- 
tween the shear viscosity rj and the jet quenching param- 
eter q for a "weakly coupled", i.e. quasi-particle dom- 
inated quark-gluon plasma. The relation associates a 
small ratio of shear viscosity to entropy density to a large 
value of the jet quenching parameter. 

The fact that 77/s saturates in the limit of strong cou- 
pling of the SYM theory but q continues to increase, sug- 
gests that the ratio T^/q may serve as a more broadly 
applicable measure of the coupling strength of a quark- 
gluon plasma. We thus conjecture that the following re- 
lations hold generally: 



1.25- 



T3 



for weak coupling. 



(25) 



^ \ q 1 for strong coupling. 

An unambiguous determination of both sides of (|25p from 
experimental data would thus permit a model indepen- 
dent, quantitative assessment of the strongly coupled na- 
ture of the quark-gluon plasma produced in heavy ion 
collisions [3]] |. 
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